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A large-S approach to the description of chiral phases in frustrated antiferromagnetic 
spin-S chains with easy-plane anisotropy is proposed. The approach is valid in the vicinity 
of the classical Lifshitz point j = \, where j is the relative strength of the frustrating 
next-nearest-neighbor coupling. The effective theory turns out to be the nonlinear <r-model 
with the additional fourth-order term in spatial derivatives. The topological term persists 
for j > j, which leads to disappearance of the chiral gapped phase for half- integer S, in 
accordance with the recent numerical results. 

§1. Introduction 

Chirally ordered phases in frustrated quantum spin chains have attracted con- 
siderable attention recently 1 ^" 9 - 1 . Nersesyan et al. 1 ' have studied the antiferromag- 
netic S = \ chain with easy-plane anisotropy and frustrating next-nearest-neighbor 
(NNN) coupling, described by the Hamiltonian 

H = ^{(S n S n+1 ) A + j(S n S n+2 ) A , {S 1 S 2 )a = S?S% + SfSl + ASfS z 2 , (1-1) 

n 

where S n denotes the spin operator at the n-th site, j > is the relative strength 
of the NNN coupling, and < A < 1 is the dipolar easy-plane anisotropy. In 
the limit j > 1, performing the Abelian bosonization and subsequently using RG 
and mean-field arguments, they have predicted the existence of a new gapless phase 
with a broken parity, which is characterized by the nonzero value of the vector 
chirality K n = ((S n x S n+ i)). This type of ordering does not break the U(l) in- 
plane rotation symmetry, it only breaks the discrete parity symmetry and thus is 
in principle perfectly allowed in one dimension, the idea probably first realized by 
Chubukov. 10 ) Except having the chiral order, this phase is characterized by the 
power-law decaying incommensurate in-plane spin correlations of the form (Sq S~) oc 
n~ ri e l ® n , where Q is very close to it in the limit j 3> 1, and n = j for S = ^. ^ 

Early attempts 2 -*' 4 -* to find this chiral gapless phase in numerical calculations 
for S = | were unsuccessful. At the same time, to much of surprise, DMRG studies 
for frustrated 5* = 1 chain 2 )' 5 ) have shown the presence of two different types of 
chiral phases, gapped and gapless. Later, chiral phases were numerically found for 
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S = ^, 7 )' 8 ) as well as for S = § and S = 2. 8 ) 

For general S, the appearance of two chiral phases was explained with the help of 
the large-,!? mapping to a helimagnet, and the qualitative form of the phase diagram 
for large S was given. 6 ) However, the way of mapping to a helimagnet used in Ref. 6 ) 
neglects the presence of the topological term and thus is in fact valid only for integer 
S, when the topological term is ineffective. 

Another theoretical approach using bosonization 9 ) suggests that the phase dia- 
gram for integer and half-integer S should be very similar, with the only difference 
that the Haldane phase gets replaced by the dimerized phase in the case of half- 
integer S. This is, however, in contradiction with the recent numerical results 8 ) 
indicating that the chiral gapped phase is absent for half-integer S. 

In the present paper I show that the difference between integer and half-integer 
S, well known for small j, 12 ) persists also for the chiral region of the phase dia- 
gram, which leads to disappearance of the chiral gapped phase and also changes 
dramatically the mechanisms of destabilization of the chiral gapless phase. 



§2. Modified nonlinear sigma model 



Consider the model (IT) for general (large) spin S. One can pass to the spin 
coherent states in a usual manner, as described, e.g., in Ref. u \ The Berry phase for 
a single spin at the site i can be expressed through the unit vector Hi parametrizing 
the coherent state: <Pi = S f dt 8t "^'^ e '- > , where e is an arbitrary unit vector. If we 
choose e*i = — H2 and e*2 = — Hi, the sum of Berry phases for two neighboring spins 
at sites 1 and 2 can be written easily in a compact form : = S n J • dt(n2 — 
Hi). Further, we select the uniform and staggered components of the magnetization, 
putting Hi = rhi + (— l) 1 ^. When we are in the vicinity of the classical Lifshitz 
point j = \, both uniform and staggered magnetization vary slowly in space, so that 
we can pass to the continuum approximation, assuming as usual that m < I. The 
effective action is readily obtained in the following form: 

A = 2t:SQ + S J dxj dtm{lx d t l) (2T) 

- S 2 J dxj - A)l\ + 2f a ml - 2eh a (d x l a ) 2 + \h a {d 2 xx l a ) 2 } , 

where Q = ^ / dt J dxl ■ (d x l x d t l) is the topological charge, and we use the notation 

e = j-l/4, f x>y = h x , y = l, f 2 = (3 + 5A)/8, h z = A. (2-2) 

For the sake of clarity, we have set the Planck constant and the lattice constant 
to unity. Making Taylor expansions of the fields, one has to take into account 
the derivatives of I up to the fourth order, since, as one can see from (2-1), the 
contribution of the second order comes with the prefactor e which becomes negative 
in the region j > \ we are interested in. Note also that we do not assume A to 
be close to 1 and thus have to keep terms like (1 — A) multiplied by m 2 or (d x l z ) 2 
etc. The uniform part m can be integrated out, and after passing to the imaginary 
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"time" r = 2iSt one obtains the following effective Euclidean action, valid in the 
vicinity of the Lifshitz point: 

AE = Jg~ J dX J dT {yV X ~ 4eK(d x l a ) 2 + - A K (d 2 xx l a ) 2 

+ 1(1- A){l z ) 2 } + i2KSQ, (2-3) 

where the bare coupling constant go = 2/5. It is easy to see that for A — > 1 and 
j — > the action (2-3) gives the well-known expression for the isotropic Heisenberg 
chain (the fourth-order derivatives become irrelevant in this limit). 

§3. Mapping to a helimagnet 

Now I will show that the action (2-3) can be further mapped to a helimagnet, 
giving the results very similar to those obtained with the help of the ansatz of Ref. 6 ) . 
Passing to the angular variables l x + il y = sinOe 1 ^, l z = cos 9, one may notice that 
for A ^ 1 the field 9 is massive, and thus can be integrated out. Putting 9 = | + ?? 
and expanding in one can rewrite the action as 



A E = A[&] + 



2<?o 



dx / dr \-{d T ^) 2 [l + (f z -2){) 2 } 



fz 

le( dx ip) 2 + \{d x vf + \{d 2 xx <p) 2 



+ i2nSQ, 



(3-1) 



m = ^-J d *J dr^) 2 - AeA{d x $) 2 + \A{d 2 x J) 2 - |(1 - A)# 

Further, introducing instead of •& a new variable £ by putting •& = ^(£e tXx + 
£*e~ lXx ), where £ is assumed to be small and smoothly varying, one can obtain the 
action A[£] where now only quadratic terms in £ have to be kept. In order to kill 
terms of the type £d x £* — S,*9 x £,, it is necessary to set A = Ao, where Ao = y/8e is the 
classical pitch of the helix in the limit The action now takes the form 



Ae = j dx J 



d,T 



±(d Tip ) 2 [i + hi + zM\ 2 

2g L 2 



+ 



1 

2m 



l 



l£l 2 V[p] 



jjdxj dr||^| 2 + \d x i\ 2 + m 2 \t\ 2 } + i2nSQ, 



+ 



where the following notation is used: 



(3-2) 



g = g (l-Z), Z=-(l-A), 



r = T\/8eZ, V[<p] 



1 



{d xl p) 2 - 8e 



3(1-^) 

2 + \(oL^) 2 - 



E 




(3-3) 



Applying the standard Polyakov-type RG, one obtains finally the effective action of 
the planar helimagnet, which depends only on in-plane angle ip: 

A[<p] = ^ J dx J dy{(d y ^) 2 + V[<p]} + i2nSQ, T cS = y/g?j, y = T^g~fg 
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9o 



9o 



1 - 



9o^ 
8irE 



ln(l+4j/m§) 



1 - 



90^ 1 + Z 
8ttE 1 - Z 



Hl + A 2 /m 2 ) 



(3-4) 



here Aq = ir is the momentum cutoff on the lattice. The problem is thus mapped to 
the classical helimagnet in two dimensions at the finite temperature T e ff . 

One may notice that we have kept the topological term in (3-4), although for- 
mally we don't have the right to do that after mapping to the planar model. Keeping 
this term is, however, important for the following discussion, if one considers singular 
configurations like vortices: In the vortex core the planar mapping becomes invalid, 
since the angle 6 deviates there strongly from ir/2. 



§4. Phase transitions 

The action (3-4) is very similar to that derived in Ref. 6 ) (although slightly 
different in detail) , and we will see that it yields exactly the same results concerning 
the transition lines. 

The transition from gapless chiral 
to gapped chiral phase is determined by 
the unbinding of vortices existing on the 
background of a state with a certain chi- 
rality. The critical "temperature" of this 
transition can be obtained by rewriting 
the action (3-4) in terms of deviation 
from the ground state with a certain chi- 
rality, i.e., setting if = Xqx + <fi. Then 
one may neglect the higher-order terms 
in derivatives of <\> and get the classical 
XY action of the form 



0.4 



0.3 



0.2 



0.1 




A = 



1 



2T c g 



J Jdxdy{(d y cf>) 2 + \ 2 (d x( f>) 2 }, 



from which the critical T e g is given by 
the equation 



T« T = (vr/2)A 



(4-1) 



Fig. 1. Behavior of the critical exponent rj of 
the in-plane correlations in the chiral gap- 
less phase (see Eq. (4-4)) as a function of 
j. Solid and dashed lines are the results 
of present approach for S = 1 and S = 2, 
respectively, and circles and squares are 
numerical results for S = 1 and S = 2 
taken from Ref. 8 ' . 



For e — > and finite (not very small) 
anisotropy I — A one has g ~ go , g ~ go , 
and (4T) translates into 

e KT = (3 + 5 A)/ Air 2 S 2 (4-2) 

exactly the result given in Ref. 6 ). In the other limit, when 1 — A —> and e is small 
but finite, one may neglect the difference between g and g, which yields the following 
equation for the transition line: 



(1 - A) KT = (16/3M5 exp{-2(7rSV2e - 2)}, 



(4-3) 
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again the result coinciding with that previously obtained. 6 ) 

The in-plane correlation function in the KT (gapless chiral) phase can be readily 
obtained and has the power-law form 



{S+(x)S-(0)) ocx 



7] = T efr /(27rA ) 



'3 + 5Z\ 



8ttS 



• 0. 
(4-4) 



The critical exponent r? increases when one approaches the transition; this behavior 
of rj is in qualitative agreement with the numerical results 8 ), see Fig. 1. One should 
mention, however, that the estimated numerical values of r\ are not universal at the 
transition boundary, while the picture advocated here would imply the universal 
KT value of r] = ^. Further studies are necessary to clarify this point. The large- 
j behavior of r], which is not accessible by the present approach, was studied by 
Lecheminant et al. within the bosonization framework, 9 ) and the predicted value 
rj(j — > oo) = 1/8S is in a good agreement with the numerical data. 8 ) 

The temperature of the transition between chiral and non-chiral phases can be 
estimated in a "solid-on-solid" approximation 13 ) by considering the fluctuations of a 
chiral domain wall. 6 ) Any such fluctuation requires formation of topological defects 
with nonzero vorticity, as shown in Fig. 2. Such field configurations, which we will 
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Fig. 2. A schematic view of the elementary "fluctuation" in a chiral domain wall: an instanton con- 
necting two domain wall states separated by the minimum possible distance n/Xo (displacement 
by one "quantum"). The instanton necessarily includes a vortex- like topological defect and thus 
has the topological charge Q = ±i. The dashed line indicates the position of the domain wall. 

call "bound vortices" , are instantons connecting states with different positions of the 
domain wall in space, and the corresponding change in the position is quantized in 
multiples of vr/Ao- The free energy a of the domain wall per unit length in the y 
direction can be written as 



a = E DW j- InZ, 



Z 



2^ dy N --- 

N=0 J Vn-i " 



dy 



N 



Ebv 



t e ff 



(4-5) 
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where L is the length of the domain wall, E^ v is the energy of the bound vortex, 
Edw is the energy of the domain wall (per unit length), No = L/do, where do is 
the minimal possible distance in y direction between successive fluctuations, and the 
positions yi of fluctuations are subject to the constraint \y\ — yi+i\ < do. The multiple 
integral in (4-5) is easily calculated to be equal to (L — Ndo) N /N\ by making the 
substitution Q = yi — (i — l)do, and after applying the Stirling formula one finally 
gets 

a = E DW - (T/do) In |l + d (cotanh(£ te /2T cfr ) - l) j (4-6) 

We cannot calculate the energy of 
the bound vortex E^ v and the scale do 
analytically. Assuming that bound and 
free vortex solutions have similar struc- 
ture apart from the core region, and cut- 
ting the free vortex solution at the dis- 
tance tt/Xq, one may estimate E^ v ~ 
7rAo ln(-7r/Ao) and do = tt/Aq- Then the 
equation a = can be solved numer- 
ically, and the corresponding solution 
T e 7 ff as a function of j for j — ► 1/4 is 
shown in Fig. 3 together with the cor- 
responding critical temperature of the 
Kosterlitz-Thouless transition. One can 
see that the two temperatures are very 
close to each other, but T 1 is slightly 
higher than T KT . 

Thus, for integer S we obtain the 
same phase diagram as suggested in Ref. 
6 ). Now, the question is what changes 
in this picture if S becomes half-integer? 

As it was mentioned before, the planar description becomes invalid near the 
vortex core, where it becomes energetically favorable to lift the spins from the easy 
plane in order to reduce the contribution of the gradient terms in the action. The 
vector I in the center of a vortex is in fact perpendicular to the plane, and thus the 
vortices present in the model are in fact not Z vortices of the purely XY system, 
but rather Z2 vortices of the Heisenberg model. The topological charge of a vortex is 
Q = —pv/2 = ±i, where p = ±1 is the sign of cos 9 in the vortex center, and v = ±1 
is the usual XY vorticity (the angle (p changes by 1-nv when one goes around the 
center). Therefore, for half-integer S every vortex obtains the effective phase factor 
e i27rSQ _ e ±«7r/2^ an j a f^ er summation over possible values of Q the contribution 
of vortices vanishes. This has the usual consequence of suppressing the Kosterlitz- 
Thouless transition, which disappears for half-integer S also for small j, together 
with the corresponding gapped (Haldane) phase. 12 -* Thus one may expect that only 
the chiral gapless phase should survive for half-integer S, in accordance with the 
numerical results. 8 ) 
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Fig. 3. The effective critical "temperatures" 
T c ff of the Ising (chiral — non-chiral) and 
Kosterlitz-Thouless (chiral gapless — chi- 
ral gapped) transitions, according to Eqs. 
(4-1) and (4-6). One can see that the re- 
gion of existence of the chiral gapped phase 
is very narrow. 
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Another, somewhat unexpected, ef- 
fect of the topological term is that it 
suppresses the fluctuations of the chiral 
domain walls as well. Indeed, since ev- 
ery elementary fluctuation (a "jump" by 
7r/Ao in space) contains a "bound vor- 
tex", the contribution of configurations 
with rtj ^ in (4-5) should vanish. That 
means that the mechanisms of destabi- 
lizing the chiral gapless phase for integer 
and half-integer S must be very differ- 
ent. At present, I am not able to suggest 
any efficient mechanism for destroying 
the chiral phase for half-integer S. One 
may speculate that high-energy config- 
urations with discontinuities may play 
some role (an example of such config- 
uration is shown in Fig. 4). Another 
possibility would be that due to interac- 
tion between bound vortices the energy 
of configurations describing the domain 
wall fluctuations with even raj will be the 
lowest when all "bound vortices" have the same topological charge Q = +5 or — | 
(note that a uniform sequence of Q means an alternating sequence of p = signcos#, 
since the vorticities v must alternate). In any case, such explanations would mean 
a strong increase of the Ising transition temperature in comparison to the integer S 
case. 
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Fig. 4. An example of configuration with dis- 
continuity which may contribute to the 
destabilization of the chiral gapless phase 
for half-integer S, see the text. The dashed 
line shows the position of the chiral do- 
main wall, and the discontinuity is indi- 
cated by the solid line. The discontinuity 
is a chiral domain wall with an additional 
7r phase jump; at both ends of the disconti- 
nuity there are "half- vortices" which carry 
the topological charge Q = ±|. 



